Skip to content

Conversation

jjokella
Copy link
Contributor

@jjokella jjokella commented Nov 5, 2024

Former dim_obs_p in add_obs_error_pdaf was actually the size of the full observation vector and is therefore better called dim_obs.

The suffix _p is used in TSMP-PDAF for PE-local variables.

Before, `clm_obserr(obs_nc2pdaf(i))` was used, where `i` is looped
over `1, dim_obs_p`. However, `obs_nc2pdaf` will then for each PE
access the first dimensions from the first PE (or more depending on
`dim_obs_p`).

The local versions `clm_obserr_p` and `pressure_obserr_p` should
contain the correct observation errors.
@jjokella jjokella marked this pull request as draft November 5, 2024 15:56
@jjokella jjokella changed the title PDAF: use PE-local observation error arrays in add_obs_error_pdaf (LEnKF) PDAF: change dim_obs_p to dim_obs in add_obs_error_pdaf (LEnKF) Nov 6, 2024
@jjokella jjokella marked this pull request as ready for review November 6, 2024 06:38
@jjokella jjokella merged commit 0450e32 into master Nov 6, 2024
1 check passed
@jjokella jjokella deleted the dev-add-obserr branch November 6, 2024 06:39
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant